Scale invariant growth processes in expanding space 
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Many growth processes lead to intriguing stochastic patterns and complex fractal structures which 
exhibit local scale invariance properties. Such structures can often be described effectively by space- 
time trajectories of interacting particles, and their large scale behaviour depends on the overall 
growth geometry. We establish an exact relation between statistical properties of structures in 
uniformly expanding and fixed geometries, which preserves the local scale invariance and is inde- 
pendent of other properties such as the dimensionality. This relation generalizes standard conformal 
transformations as the natural symmetry of self-afhne growth processes. We illustrate our main re- 
sult numerically for various structures of coalescing Levy flights and fractional Brownian motions, 
including also branching and finite particle sizes. One of the main benefits of this new approach 
is a full, explicit description of the asymptotic statistics in expanding domains, which are often 
non-trivial and random due to amplification of initial fluctuations. 
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Scale invariant structures resulting from fractal grovifth 
processes are abundant across nature [l], Q ■ Examples 
include diffusion- limited aggregation 0], river basins 3, 
and self-affine domain boundaries forming behind grow- 
ing fronts for spatial competition models Q. While 
first results appeared already 30 years ago, the field 
continues to be of interest [7| with recent applications 
in microbial growth @, Q. In many cases these struc- 
tures can be modelled as trajectories of locally inter- 
acting particles — a picture that wc adopt in this Letter. 
The overall geometry has a strong impact on growth pro- 
cesses. A dramatic example is viscous fingering, where 
in constant width channel geometry a stable Saffman- 
Taylor finger of fixed shape propagates [l^l , while in ra- 
dial geometry a continuously tip splitting branched struc- 
ture emerges ll|, |l2| . In biological growth spatial range 
expansion is often coupled to drift and competition in 
the genetic pool |13j . and is rccongnizcd to have major 
influence on the gene pool of natural populations [l4| . 

In this Letter we show how the effect of the overall 
geometry in many directed growth processes can be cap- 
tured elegantly in terms of a time dependent metric. We 
view growing domain boundaries as space-time trajec- 
tories of particles moving on the growth front, which is 
expanding in many interesting cases. A natural example 
within the scope of this Letter is isoradial growth in two 
dimensions, such as domain boundaries of competing mi- 
crobial species in a Petri dish 0] . While cosmology is an 
obvious example, there has been recent interest in non- 
constant metric also in thin sheets 15-l3]- Our results 
are applicable to the formation of stochastic patterns and 
structures in a very general setting, including diffusion 
processes with time-dependent diffusion rate (i.e., tem- 
perature) [l8-21|, in cosmologically expanding space p^ . 
or on a biologically growing substrate. 



In particular, we consider self-affine space-time tra- 
jectories of particles under spatially homogeneous but 
time dependent metric, and map those into more eas- 
ily tractable systems with constant metric. The map- 
ping depends only on the local scale invariance exponent 
of the trajectories, and works directly for local interac- 
tions which do not involve a length scale, such as an- 
nihilation or coagulation of point particles. Branching 
or exclusion/reflection of finite size particles can also be 
treated after mapping the interaction length scales ap- 
propriately. This provides a natural extension of con- 
formal maps to generalized self-affine growth processes, 
and we show how this leads to an exact description of 
the non-trivial asymptotic statistics of growth structures 
in expanding domains, which is one particularly striking 
consequence of this new approach. 

To describe our results in the most illustrative setting, 
we consider the growth of self-affine structures (e.g., do- 
main boundaries) in isoradial geometry in two dimen- 
sions. These structures consist of directed "arms" , which 
can be interpreted as locally scale invariant space-time 
trajectories of point particles moving in an expanding one 
dimensional space with periodic boundary conditions. 

Consider an isotropic radial structure growing from an 
initial disk with radius tq, shown in Fig.[lja) for an exam- 
ple of radial coalescing Brownian motions, where also the 
following notation is illustrated. We describe each arm 
by the displacement along the perimeter of the growing 
circle 



Yr e [0, 27rr) with r > tq 



(1) 



as a function of the radial distance r; directed radial 
growth means that this is possible. In the increment 



dYr = Yr dr/r + dYr 



(2) 
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FIG. 1: (color online). Expanding radial growth structure and 
the same structure on a fixed domain with periodic boundary 
conditions, illustrated for the case of coalescing random walks 
(7 = 1/2). The notation for the displacements Yr ^ and Xh 
^ is illustrated in blue. The distribution of the rescaled ra- 
dial structure at radius r is identical to the distribution of the 
fixed domain structure at height h{r) as given by (jS}, indi- 
cated by dashed red lines. This mapping (plotted in Fig. [2]) 
has a finite limit h.y{oo) for 7 < 1, indicated by a dashed black 
line. Parameters are L = 100 with ro = L/2tt, unit diffusion 
coefficient and initially 100 arms. 




FIG. 2: (color online). The mapping h{r) between an expand- 
ing radial growth structure and the same structure on a fixed 
domain as given in Q in units of ro, see (lllf) . Due to local 
equivalence of the two processes, h{r) w r — ro for r ~ ro. 
The different geometries affect the behaviour at large r, in 
particular h{r) has a finite limit h^{oo) for 7 < 1 (|10|l and 
diverges for 7 > 1. The asymptotic behaviour is indicated by 
dotted lines, except for 7 — 2/3 which is off the figure. 



the first term is due to the stretching of space, and the 
second corresponds to the inherent fluctuations encoding 
the local scale invariance of the arms. Instead of ra- 
dial coordinates (F, r), the arms can also be represented 
in modified polar coordinates {X,h): the polar angle is 
multiplied by vq and denoted by X^, which is in a fixed 
periodic domain: 



Xh e [0, L) with h>0, 



(3) 



and the relation between r and h will be determined 
shortly. The choice L ~ 27rro enables matching the initial 
conditions between Xh=o ^md Yr^ro- This implies 



which using Eq. yields for the increments 
dXh =—dYr. 



(4) 



(5) 



We impose that the mapping between expanding and 
fixed geometry preserves the relevant local structure of 
the object (analogously to conformal invariance), which 
in our case is given by local scale invariance of the arms 



dXh - (dhy and dYr - (dr) 



(6) 



with 7 > 0. For example, diffusive fluctuations corre- 
spond to 7 = 1/2 (see 123| for related results), and for 
ballistic displacements of the arms 7 = 1. Other values 
are related to sub- or superdiffusive behaviour, such as 
7 = 2/3 for domain boundaries driven by a surface in the 
KPZ universality class [i,[2i,[25|. 



This leads to a relationship between h>0 and r > ro 
via 



dh _ UXh\^'^ fra^^h 
dr ~ \ dYr 



(7) 



where multiplicative prefactors, which are equal in ^ 
although not indicated, drop out. Integrating yields 



h{r) = I 

for all r > rg 
tial condition Y, 

r^y di|t. ^ 



ro hi(r/ro) 



7 = 1 



(8) 



h{r} 



For a single arm the matching ini- 
= Xq leads to identical distributions 

for all r > ro. Our main result is 



now that the same holds for the entire growth struc- 
tures which are characterized as collections of arms {Yr} 
and {Xfi}, with the independent variables linked through 
h = /i(r): 



forallr>ro 



(9) 



provided that the arms interact only locally. Examples 
of such interactions include coagulation or annihilation, 
and wc discuss how this can be generalized in more de- 
tail below. Figure [1] illustrates this correspondence for 
coalescing Brownian trajectories. 

Properties of the mapping. To leading order h{r) sa 
r — ro for r close to ro, since locally the fixed domain 
and the radial models are equivalent. The effect of the 
different geometries enters in the non-linear behaviour of 
h(r) for larger values of r, in particular for < 7 < 1 we 
have 



h^{oo) = lim h{r) 



1 



1-7 



ro < 00 . 



(10) 
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FIG. 3: (color online). Numerical demonstration of the mapping ((8]) between radial geometry (o on all panels) and fixed 
domain (x): (a), (b) number of surviving arms (A'') and their root mean square distance (-D^) (|15|l for Levy flights H13|) with 
7 = max{l/a, 1/2}; (c), (d) same observables for fractional Brownian motion (|14p : (e) number of surviving arms for Brownian 
motion (7 = 1/2) with finite particle size d, + symbols indicate uncorrected, x corrected particle radii (see text); (f) radial 
density of arms {N)/r for branching-coalescing Brownian motion with fixed branching rate Rr (|16p in the radial geometry. 
For panels (a)-(e) the horizontal axis is h or h{r), with h{oo) indicated, while for panel (f) the horizonal axis is r or r{h). The 
asymptotic scaling laws [black dashed lines on panels (a)-(d)] break down when (A*") ^ 1. Panels (a)-(d) and (f) are in 1 + 1 
dimensions, initially 100 arms, L = 100 with ro = i/27r, while (e) is in 2 + 1 dimensions with ro — 20. 



This observation is particularly interesting for coalescing 
or annihilating structures, which exhibit an absorbing 
state in a fixed geometry with one or no arms remaining 
as /i — 00. Such structures often occur in neutral mod- 
els for competition in spatial populations Q , and the 
absorbing state corresponds to fixation of the model in 
one of the initial types. By standard arguments the time 
to fixation scales as L^/'' ~ fg^^'i which is much larger 
than h{oo) ^ rg for large systems. For 7 < 1 we not only 
confirm the previous (intuitive) result that there is no fix- 
ation in expanding populations in the neutral case, but 
also give explicitly the spatial distribution of the surviv- 
ing types at large radii r ^ 00 as {Xfi(^ao)}- The process 
Xh is much easier to simulate than Yr , and in many cases 
there also exist theoretical predictions [iij . 

In Fig. [2] we plot the mapping in convenient dimen- 
sionless variables r' = r/ro and h' = h/rQ, so that 



/i'(r') = J 1-^ ' ^ 



- (1/r')" 
log(r') 



,7^1 
,7-1 



(11) 



for all r' > 1. For 7 = 1 we recover the generic confor- 
mal map from the exterior of the unit circle to a strip, 
whereas for 7 7^ 1 the mapping provides a natural gener- 



alization to self-affine processes. Note that for general 7, 
h'{r') ~ log^(r') is the g- logarithm with a = I/7 known 
from non-extensive statistical mechanics j27l |. which can 
therefore also be interpreted as a generalization of con- 
formal transformations. 

It is instructive to consider the mapping also for inward 
growing radial structures, where r < rg (i.e., r' < 1), 
which formally leads to negative heights h < 0, corre- 
sponding to a fixed domain structure growing downward. 
Observing the general relation 



h'{l/r') = -{r')^ h'{r') for all 7 > , 



(12) 



all phenomena for such structures can be entirely un- 
derstood by studying outward growing ones. Note that 
in contrast to the expanding case now all sub-ballistic 
structures lead to fixation since |/i'(?'')l c)o as r' — 0, 
whereas super-ballistic structures will have a non-trivial 
limit. First results on inward growing radial structures 
have been obtained in [28| and our approach provides a 
framework for a better understading of those which is 
explained in detail in future work [26| . 

Validity and locality. The mapping is based purely on 
a conservation of local scale invariance of the structure. 
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Therefore it is not surprising, that the mapping can be 
shown to hold rigorously for processes which are fully de- 
termined by their local structure, namely processes with 
independent increments such as Brownian motion and 
self-similar Levy processes [l^ . On the other hand, there 
are other self-similar processes with the same local scale 
invariance but more complicated temporal correlations, 
such as fractional Brownian motion (fBm) [sO] • The cor- 
relations will influence the mapping and it does not hold 
in general for such processes. Using fractional stochas- 
tic calculus, one can derive a similar mapping for the 
particular model of fBm, which leads to a more complex 
expression which is numerically very close to ([5]). This 
derivation is beyond the scope of this letter and is dis- 
cussed in detail in [2^ . 

In Fig. [3] we illustrate the validity of the mapping for 
self similar Levy flights, which are defined via indepen- 
dent stationary increments with an a-stable jump size 
distribution 

F{Xh+ - X,, = x) ^ C\x\-^^+°'^ (a>0), (13) 

as well as fBm, which can be characterized as a Gaussian 
process with covariances 



(14) 



Levy flights have local scale invariance parameter 7 = 
max{l/a, 1/2}. They are super-diffusive and have non- 
continuous paths for a < 2, and scale diffusively for a > 2 
where the jump size has finite variance. fBm can be 
super- or sub-diffusive and is not Markovian, but still 
the mapping (jS]) works very well also in that case. In 
Fig. |3] we compare two statistics for coalescing interac- 
tion: the average number of arms {N), and the total 
mean squared distance between neighboring arms, as a 
measure for their spatial distribution. For fixed geometry 



N{h) 

E 



X 



(*+l) 



-X 



(15) 



with an analogous Dpi{r) for radial geometry. Plotting 
the fixed and circular data against h and h{r), respec- 
tively, we obtain a data collapse. The power-law pre- 
dictions for the fixed system in panels (a) to (d) can be 
derived easily by standard mean- field arguments [3l|, [s^l . 

A natural step to include non-local interactions is to 
introduce a particle size. For simplicity we consider 
isotropic shapes with diameter d, i.e., particles coagu- 
late or annihilate already at a non-zero distance d. As 
long as the diameter is much smaller than the macro- 
scopic length scale in the system, d ^ tq, the corrections 
introduced are small. Still they can be taken into ac- 
count exactly by comparing the radial system with the 
fixed domain one, where the particle diameter decreases 
as T^*^- In Fig.[3l^e) we show both cases, with and with- 
out this correction, for coalescing Brownian motions. We 



see that for small d the mapping still works very well even 
without corrections. Unlike all other numerical data pre- 
sented in this paper, this one is for an expanding sphere 
in 2 -|- 1 dimensions. Finite range interactions are partic- 
ularly important in higher dimensions, where coalescence 
or annihilation of point particles does not strictly occur, 
they only get arbitrarily close to each other. The map- 
ping is independent of the dimension, as discussed below. 

Another natural interaction included in growing struc- 
tures is branching. This is not a purely geometric in- 
teraction but has its own characteristic rate R, which 
introduces a time scale in the system. For the mapped 
processes to have the same statistics we require that the 
number of branching events Ap{dh) in the fixed domain 
model during a time interval dh is the same as Aii{dr) for 
the corresponding radial system. This implies a relation 
between the branching rates 



Rr _ Afl(dr)/dr _ dh 
'r^ ^ AF{dh)/dh ~ rf^ 



(16) 



which is {tq/tY/^. Thus to understand the density of 
branches Nr (r) /r in a radially growing system with fixed 
branching rate Rr, one has to compare to a fixed do- 
main system with increasing branching rate Rpih) ^ 
RurihY /tq, where r{h) is the inverse of (O. Note that 
this rate diverges asr— >ooor/i— >/i(oo). The density 
of branches for three different branching rates is shown 
in Fig.[2)Jf) for Brownian motions with 7 = 1/2. 

Generalized geometries. Our results can be directly 
generalized to an arbitrary time dependent domain of 
size L(t) with homogeneous metric. We obtain 



m = f 

Jo 



L(0)^^l/7 



ds 



(17) 



analogously to Eq. ([5]) . For example one can study expo- 
nentially increasing domains, which is analogous to struc- 
tures with exponentially decreasing diffusivity. These 
have been studied in detail for single random walks 
21 1 and are used in simulated annealing (33| . 

In n -|- 1 dimensions, where n is the spatial dimension- 
ality, our method applies directly if the scale invariance 
holds in all spatial directions i ~ 1, . . . ,n 



dX, - (dh)'' and dY, - (dr) ' 



(18) 



It is possible to have anisotropy (possible i-dependence 
of the multiplicative factors which are not indicated), but 
7 should be identical in all directions. Then the mapping 
P?)) stays exactly the same. 

Summary. We have demonstrated that a large class 
of locally scale invariant, directed complex structures 
growing in radial or general increasing geometries can 
be mapped to structures in fixed domains, which are 
simpler and for which exact results are often available. 
This approach provides an elegant and remarkably simple 
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way to understand various phenomena related to time- 
dependent metric, such as the effect of range expansions 
in competitive biological growth. A particularly strik- 
ing example is a full description of the limiting statistics 
of radial competition interfaces. Further examples and 
technical aspects are discussed in more detail in jiij . 
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